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ABSTRACT: We describe a technique for constructing the effective chiral theory for 
quenched QCD. The effective theory which results is a lagrangian one, with a graded sym- 
metry group which mixes Goldstone bosons and fermions, and with a definite (though slightly 
peculiar) set of Feynman rules. The straightforward application of these rules gives auto- 
matic cancellation of diagrams which would arise from virtual quark loops. The techniques 
are used to calculate chiral logarithms in fx/ fir, rrhr, f^K-, an d the ratio of (ss) to (uu). 
The leading finite-volume corrections to these quantities are also computed. Problems for 
future study are described. 



I. Introduction and Motivation 



The quenched approximation to QCD, in which virtual quark loops are neglected, is a 
necessary evil in lattice QCD simulations and will be with us for the foreseeable future. Even 
with the proposed QCD Teraflop Machine, 2 the quenched approximation will be needed to 
approach the crucial corners of the parameter space: large volumes, physical quark masses, 
and the continuum limit. We therefore need to learn as much as possible analytically about 
the quenched approximation in order to have good control over the systematics of such 
calculations. 

In the full theory, chiral perturbation theory (ChPT) is a key analytic tool. It gives: 

• The detailed form of the approach to the chiral limit. The universal terms ("chiral 
logarithms") can be calculated order by order in the loop expansion. Comparison with 
this expected chiral behavior provides, for example, a crucial check of lattice weak matrix 
element calculations. 

• The leading finite-volume corrections at large volume. 3 As the lightest particles, the 
pseudoscalar mesons clearly control these corrections; ChPT is simply the effective theory 
of their interactions. 

It is therefore clear why one would like to have a ChPT corresponding to the quenched 
approximation. In fact, there have been several previous attempts to calculate quenched 
chiral logarithms. Morel 4 and Sharpe^ use the strong coupling and 1/d expansions; Kilcup 
et a/. 6 and Sharpe 7 use the quark-flow approach (see below). The papers by Sharpe in 
particular emphasize the importance of quenched ChPT and discuss several of the key issues 
(in particular, the problems caused by the r] ). A preliminary version of the current work 
has been presented in ref. 8. 
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II. Quark-flow Approach 



In this approach, one starts with ordinary ChPT for full QCD and writes down all 
meson diagrams which contribute to the process of interest. To each meson diagram one 
then associates one or more quark-flow diagrams in QCD. Next, one eliminates all those 
quark-flow diagrams which have virtual quark loops. Finally, one attempts to reinterpret 
this elimination as conditions on the meson diagrams. Note that, in the case where more than 
one quark-flow diagram corresponds to a given meson diagram, it is by no means obvious 
that the final step can always be performed. We have been able to carry it through, more- 
or-less satisfactorily, in simple cases (see below), but have been unable to prove — within 
the context of this approach — that it can always be done. 

In order to go back and forth between quark flow diagrams and meson diagrams, the 
natural basis to use is the qq basis. In the neutral sector, this means that one works with uu, 
dd, and ss states rather than 7r°, rj, and rj' . The latter basis is convenient in the full theory 
since one can treat the rf mass as "large," decouple it, and work only with 7r° and rj. This 
turns out not to be possible in the quenched theory. In full QCD the rf gets the singlet part 
of its mass (= /z) through the iteration of quark loop diagrams joined by gluons (see fig. 1). 



Fig. 1. Quark flow diagrams for the rj' propagator in full QCD. 

In the approximation where the rf mass is much greater than the masses of the octet 
mesons, the rj' decouples and may be neglected. In the quenched approximation, on the 
other hand, only the first two diagrams in fig. 1 survive, and only the second diagram (the 
"two-hairpin" diagram — fig. 2) depends on 11. 
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Fig. 2. The "two -hairpin" diagram, the only diagram which distinguishes the singlet 
from the octet meson propagator in the quenched approximation. 

The "two-hairpin vertex" in fig. 2 is ~ /A Since the vertex is not iterated, /z 2 appears 
in the numerator, not in the denominator, of the 77 propagator. Thus the 77 cannot be 
neglected in the quenched approximation. 

In many cases, it is immediately clear which full ChPT diagrams should be dropped 
in the quenched approximation. For example, consider the lowest order correction to a n + 
propagator: a meson tadpole, shown in fig. 3. 



Fig. 3. The one-loop contribution to the pion propagator for full QCD, 
with 71 + , 7T°, K + , , rj, and rf on the loop. 

When the tadpole is a K + , then the diagram must be absent in the quenched approx- 
imation, since the s quark in the K + is not present in the external states and must come 
from a virtual loop. When the tadpole is itself a 7r+, however, the situation is less clear. If 
the s quark of the previous K + tadpole is replaced by a d, then again the diagram is absent 
in the quenched approximation. But there is now a second possibility: the valence quarks, 
themselves, could make the tadpole as in fig. 4. 
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Fig. 4- A possible valence quark contribution to the pion propagator. 

The vertex in such a diagram is a meson-meson scattering vertex with no quark exchange. 
It turns out that such a vertex vanishes at 0(p 2 ) in ChPT, although we have never been 
able to prove this to our complete satisfaction within the quark flow approach.* Thus n + 
tadpoles are also absent in the quenched approximation. Indeed, the only correction to the 
quenched n + propagator at this order comes from an rf tadpole with a single two-hairpin 
vertex, fig. 5. 



Fig. 5. The quark flow diagram for the one-loop contribution 
to the pion propagator in quenched QCD. 



The end result of this approach can be described as a "lagrangian + rules." The la- 

grangian is the ordinary chiral lagrangian corresponding to full QCD. The rules give the 

For the number of flavors Np > 4, it is easy to show that this vertex vanishes. In that 
case, one can choose all 4 participating quarks to be different and thereby make a unique 
correspondence between the quark vertex and a meson vertex. Examination of the trace 
structure of the 0(p^) chiral lagrangian then immediately gives the desired result.^ However, 
the proof in this context for Np = 3 escapes us, though the vertex certainly does vanish, as 
can be seen by working backward from the known result derived in the lagrangian framework 
of section III. 
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weighting of the diagrams and prescribe how to replace r] contributions by two-hairpin di- 
agrams. However, we find this approach unsatisfactory for two reasons. First of all, it is 
difficult to make the application of the rules routine. As mentioned above, it is not obvious 
that one can always interpret the elimination of quark-flow diagrams with virtual loops as 
conditions on the mesons diagrams. One is, at the minimum, forced to prove the vanishing 
of various vertices (such as the meson-meson no-exchange vertex). Not only are convincing 
proofs elusive, but new processes may bring up new such vertices, so it is never clear that 
all the necessary proofs have been produced. 

A second, more fundamental, problem with this approach is that the presence of the 
"rules" implies that one does not have a true lagrangian theory. This means, for example, 
that the invariance of the physics under field redefinitions is not guaranteed. Such redefi- 
nitions are needed to reduce the number of terms in the lagrangian involving the 77 — see 
ref. 10. (These terms are not constrained much by symmetry because of the anomaly.) 
Similarly, the cancellation of the quartic divergences in ordinary ChPT is guaranteed by the 
chiral invariance of the lagrangian and the measure. It is not clear (at least to us) whether 
the rules automatically respect this cancellation. 

We therefore turn now to an alternative approach to quenched ChPT which gives a true 
lagrangian framework and makes the calculation of quenched chiral logarithms routine. 

III. A Lagrangian Framework 

We start with QCD. To make a lagrangian which describes the quenched approximation, 
we take the ordinary QCD lagrangian and add, for each quark q a (a = u, d, s), a scalar (ghost) 
quark q a with the same mass. 4 The ghost determinant then cancels the quark determinant. 
Of course, the resulting theory is not unitary in the quark sector; that is acceptable since 
the quenched approximation is not unitary. 
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Assuming that quark confinement still holds, the low-energy effective theory for this 
quenched QCD lagrangian may now be constructed. It will describe the interactions of 
all possible pseudoscalar bound states of quarks or scalar quarks with their antiparticles: 
ordinary qq mesons (tt, K, . . .) which we denote, generically, by 0; ghost qq mesons denoted 
by 0; and fermionic mesons qq and qq denoted by \ an d X > respectively. 

As in ordinary ChPT, the symmetries at the quark level determine the form of the inter- 
actions among the mesons. The symmetry is U(3\3)i x U(3\3)r, where U(3\3) is "almost" a 
U(6) among u, d, s, u, d, s, but has a graded structure since it mixes fermions and bosons. 11 
If we write a matrix U G U(3\3) in block form as 

u = (d £)■ « 

then A and B are 3x3 matrices of commuting numbers; C and D, of anticommuting. 
Unitarity is defined as usual: WU = I. Hermitian conjugation (f) also has the usual 
definition (complex conjugation of the usual transpose), but complex conjugation is defined 
to switch the order of anticommuting variables: (ej^)* = e^e^. There is also a cyclic 
"supertrace" defined by str(U) = tr(A) — tr(B), and a "superdeterminant," sdet(U) = 
exp strlnU, with the property sdet{U\U2) = sdet{U\)sdet{U2)- Explicitly, 

sdet{U) = det(A - C B' 1 D) / det(B) . (2) 

Now define the Hermitian field $ and the mass matrix M. by 

where 

fm u \ 

M = m d , (4) 

V m s J 

is the usual quark mass matrix. Note that, to lowest order in M, these ChPT quark masses 
are the same as those of QCD. 
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The unitary field E = exp(2«$//) transforms as E — > UjT^W^. The lagrangian invariant 
under the full U{3\3) L x C/(3|3) fi is then 

Anv = —strid^T)) + u str(7WS + ME*), (5) 
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where / and t> are as yet undetermined bare coupling constants. This looks very much like 
ordinary ChPT. 

The anomaly breaks the symmetry group down to SU{3\3)i x SU (3\3)r x £/(l). The 
anomalous field is $o = (v' ~ ^O/v^, where the minus sign comes from the relative minus 
sign between boson and fermion loops. Under the reduced group, $q ^ strlnE = InsrfetE 
is invariant, so arbitrary functions of $q can be included in the full lagrangian, C However, 
in the current framework one can redefine E to simplify C, much as in ref. 10. The result is 

C = Vi($o)str(<9 M E<^Et) + V 2 ($ )str(MX + MT$) 

-Vb($o)+V5($o)(^$ ) 2 , (6) 
£($ = 0) = Anv, (7) 

where the functions can be chosen to be real and even by making use of the freedom 
allowed by field redefinitions. In ref. 10 a different choice is made: V5 is set to but V2 
is kept complex. The potentials V3 and V4 from ref. 10 are not needed for the purposes of 
this paper and have been dropped. Note that the notation in eqns. (6) and (7) is slightly 
different from that used in ref. 8. For the purposes of this paper, we need only the quadratic 
terms in £$ . We have 

C = £ imr + a(d fl Qo) 2 -n 2 Q% + ..., 

a = \/ 5 (0), (8) 
/i 2 ee (l/2)Vb"(0). 

One can now calculate straightforwardly with C Note that because of the minus sign 
in the definition of str, some of the fields will have negative metrics. An unusual feature 
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occurs in the rj', fj' sector: terms from Vo and V5 have a different matrix structure from those 
of £i nv , and one cannot diagonalize the quadratic lagrangian in a momentum independent 
way. This leads us to treat the quadratic terms from Vq and V5 as vertices. Iterations of 
these vertices on the same line then automatically vanish due to cancellation between the 
7/ and the negative metric fj' . This is a manifestation of the fact that the iteration of the 
two-hairpin vertex is forbidden in quenched QCD. When m s 7^ m another peculiarity occurs: 
the 7T° is the only well-behaved neutral particle. The propagators of the orthogonal states 
do not have simple-pole structures. When also m u 7^ m^, even the 7r° propagator becomes 
ill-behaved. 

Because of the unusual structure of the neutral sector, it is convenient, both in the 
formalism and in actual computations, to write the neutral meson propagators in the basis 
of the states corresponding to uu, dd, ss and their ghost counterparts. As mentioned above, 
this is unlike the case of full QCD, where, due to the singlet part of the rj mass, the 
propagators in this sector are diagonal in the 7r°, 77, 7/ basis. 

Since we have a true lagrangian theory, the symmetry should guarantee that any quartic 
divergences in the diagrammatic expansion will be cancelled by contributions from the mea- 
sure, just as in the full theory. We have explicitly checked this in the S77(l|l)_£, x SU(1\1)ji 
case. The SU(3\3)i x SU(3\3)r case is considerably more complicated; however, it turns 
out that no quartic divergences appear in any of the calculations presented below. 

The vertices from Vo and V5 can in principle appear more than once in a diagram if they 
occur on different lines. However, it is our philosophy to treat the parameters /j and a as 
small. This is certainly true in the 1/N C expansion. (Recall that in full QCD the 7/ gets its 
singlet mass at order l/iV c . 12 ) Moreover, it appears that the real expansion parameters are 
a/3 and /i 2 /3 (see below). To estimate the size of the one-loop corrections, one may take 
a = 0, neglect 77-7/ mixing, and use the physical rf mass. One gets /i 2 /3 ~ (500 MeV) 2 ~ 
which leads one to expect that quenched ChPT should be roughly as good as full ChPT for 
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the kaon. 



IV. Results and Conclusions 



We have calculated, at one-loop, m^, mx, fir, fx, ( uu ), an d ( ss )- m t ne isospin limit 
{m u = rrid = m) and at infinite volume, we get: 

{m ^f = m\ (l + ^ (|A> - £ + §™» + (£ - |m?) In^/m?; 
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(^ l0 ° P ) 2 = -1,(1 + ^2 (f A 2 + - j4) 1-(A 2 /^) 



3 3 K K ) 2{m\-m\) \ m| 

/^ 100P = /, (9) 
f K \ 1 - ]oap _,, 1 2 f im K ^ml(2m K -ml) (2ml 



m ( MM )l-looP = _I( m l- lo °P)2 / 2 



^ - ^(2m| - ml^j ln (A 2 /(2m| - m 2 )) + 2L(2m 2 K - mj) - ^ 
where A is the cutoff, and 7?%, and / are the bare parameters: 



2 _ 8t;m _ Mm s + m) 
m ^r - -^2~ ' m if - j2 • ' ' 



It should be noticed that, except for the A 2 terms, a terms are actually higher order in a 
combined expansion in 1/N C and M. This implies that, apart from quadratically divergent 
terms, we may set a = in eq. (9) systematically. 

Note that in the quenched approximation the ratio fx/ fir is finite at one loop, unlike 
the full theory where this quantity contains a logarithmic divergence. In fact, if we consider 
ratios in which the quadratic divergence cancels, and then set a = as argued above, the 

10 



ratios (m^ loop /m^- loop ) 2 and ((ss) / (uu)) 1 loop are also finite. Expressed in terms of the 



bare quark masses we have 



1— loop\ 2 , 
m K \ m + m s f pi z /3 



(a») 1 - loop . " 2 /3 w , , 

1 — i = 1 o— o InlTOo/m). 

(mm) 1 - 100 ? 8tt 2 / 2 

Despite the fact that many of our one-loop results are finite, they are not quantitative 
predictions because the terms in the C(p 4 ) lagrangian (Gasser and Leutwyler's 10 L^s) may 
also contribute. In other words, we have computed the "chiral logarithms" only, not what 
are usually called the "finite terms" (which are always uncomputable in ChPT). One would 
need in general to take further ratios of physical quantities to eliminate such uncertainties. 

Using a = and u as estimated above and neglecting the "finite terms," {f k I f^) l ^°° V — 
1.07, indicating that quenched ChPT is working well. Note however that (r?% loop /m 7r ) 2 = 
1.5 for A = lGeV and ((ss) / (-u-u)) 1_loop = 0.4, although these ratios are not directly physical, 
and the large corrections should perhaps not be worrisome. 

We have also computed the leading finite volume corrections to the above results. The 
calculation is straightforward: we simply replace the infinite-volume meson propagators by 
their finite-volume counterparts. We find: 

A ((rn)- 1 """) 2 ) = 0, 
A (/i- l00p ) = 0, 



,2 



(12) 



271 e -rn„L^ 



A (m^ss) 1 - 100 ^ =0, 
where A (quantity) denotes leading order finite volume corrections to be added to the infinite 
volume one- loop expressions for the quantities given in eq. (9). L is the spatial size of the box. 
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We have assumed periodic boundary conditions* and T»i, where T is the temporal size of 
the box. In addition we have neglected terms of order e - m ^(v / 2i) an( j f or( j er e -m K L^ Thus 
the leading corrections come only from pions propagating from the closest periodic images of 
the original box. Finally we have taken only the leading contribution to the pion progagator 
and neglected terms of order (m 7r L)~^e~ TO7r ^. Note that the terms we have neglected may 
very well not be small in many current lattice simulations. In such cases, the exact one-loop 
finite volume corrections (computed by using the exact finite- volume propagators) should be 
used. 

Some comments on our results and directions for future work: 

• The absence of chiral logarithms in seen in ref. 5 is presumably a feature of the 
leading term in the 1/d expansion. Indeed, the 7/ diagrams which give such logarithms 
are mentioned in ref. 7. 

• Many of these results (e.g., Uk/U) 1 ^°° V , ((ss)/(mm)) 1 - 1oo P, or m^ loop /m^) blow up 
as m — > 0. This is an IR effect coming from the double pole of the two-hairpin diagram 
and is absent in the full theory where the vertex is iterated. It is not clear at this point 
whether this is a sickness of the quenched approximation or only of the current quenched 
chiral expansion. 

• A Gasser-Leutwyler 10 program for quenched ChPT at one-loop is possible: we expect 
there to be interesting numerical relations involving only computable (on the lattice) 
quantities. Such relations could give quantitative insight into the effects of quenching. 

• The techniques described here can be easily used to generate the effective theory for 

a QCD in which the quark loops are not neglected, but the masses of valence and 

virtual quarks are not identical. Such an effective theory is relevant to many "full 

* In order to preserve the graded symmetry, the ghost quarks must have the same boundary 
conditions (antiperiodic or periodic) as are chosen for the quarks. All mesons, including the 
fermionic ones, will thus have periodic boundary conditions. 
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QCD" simulations. Similarly, one can generate effective theories which correspond to 
the quenching of some, but not all the light quarks. 

• It seems to be straightforward to extend these ideas to the calculation of the chiral 
logarithms in weak matrix elements. 
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